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Abstract. One of the central challenges in kinetic theory is the derivation 
of macroscopic evolution equations — describing, for example, the dynamics of 
an electron gas — from the underlying fundamental microscopic laws of clas- 
sical or quantum mechanics. An iconic mathematical model in this research 
area is the Lorentz gas, which describes an ensemble of non-interacting point 
particles in an infinite array of spherical scatterers. In the case of a disor- 
dered scatterer configuration, the classical results by Gallavotti, Spohn and 
Boldrighini-Bunimovich-Sinai show that the time evolution of a macroscopic 
particle cloud is governed, in the limit of small scatterer density (Boltzmann- 
Grad limit), by the linear Boltzmann equation. In this lecture I will discuss 
the recent discovery that for a periodic configuration of scatterers the linear 
Boltzmann equation fails, and the random flight process that emerges in the 
Boltzmann-Grad limit is substantially more complicated. The key ingredient 
in the description of the limiting stochastic process is the renormalization dy- 
namics on the space of lattices, a powerful technique that has recently been 
successfully applied also to other open problems in mathematical physics, in- 
cluding KAM theory and quantum chaos. This lecture is based on joint work 
with Andreas Strombergsson, Uppsala. 



1. Introduction 

An important cornerstone in mathematical physics is the problem of deriving 
macroscopic evolution equations from first principles, i.e., the microscopic laws of 
motion governed by quantum theory, or (to simplify) Newton's laws of classical me- 
chanics. The subject has its origin in Boltzmann's revolutionary vision formulated 
more than a century ago, and it is perhaps surprising that still today there is no 
complete understanding of his most fundamental model, a dilute gas of hard spheres: 
Lanford's seminal work [29 establishes the validity of the Boltzmann equation only 
for times that are a fraction of the mean collision time. Even simpler models, such 
as a gas of non-interacting particles in a fixed array of scatterers studied first by 
Lorentz in 1905 |31j (the Lorentz gas), are difficult to analyze and lead to new and 
unexpected macroscopic phenomena. This lecture will focus on the Lorentz gas and 
report on recent joint work with Andreas Strombergsson [37], [3H], [39], [10] in the 
case of a crystal with periodic scatterer configuration. The exciting aspect of our 
findings is that the kinetic transport equations that emerge in the limit of small 
scatter size (the Boltzmann-Grad limit) are new, and distinctly different from the 
answer for disordered scatterers. 

In his original work, Boltzmann considered the case of a dilute gas of hard spheres 
with clastic collisions; dilute refers to the limit of small particle density, i.e., the ratio 
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Figure 1. The Lorentz gas for a random scatterer configuration. 

of particle radius and mean separation tends to zero. In the case of a gas of non- 
interacting particles in an array of fixed scatterers — the Lorentz gas — the dynamics 
can be reduced to a one-particle motion (Fig. [l]). We denote by q(t) £ R d the 
position and by v(t) £ M. d the velocity of our particle. The correct kinetic limit is 
then obtained by taking the scatterer radius p —> 0. A simple dimensional argument 
shows that in this limit the mean free path length should scale like p - ^ -1 ), i.e., the 
inverse of the total scattering cross section of an individual scatterer. This suggests 
to rescale the length units by introducing the macroscopic coordinates 

(1) (Q(t),V(t)) = {p d - l q{p-^- l h),v{p- {d - l) t))- 

This rescaling of length and time is commonly referred to as the Boltzmann-Grad 
scaling, and the corresponding limit p — > as the Boltzmann-Grad limit. The time 
evolution of a particle with initial data (Q, V) is then described by the billiard flow 

(2) (Q(t),V(t)) = $l(Q,V). 

Since the speed of our particle is a constant of motion we may assume without 
loss of generality that || V|| = 1. For notational reasons it is convenient to extend 
the dynamics to the inside of each scatterer trivially, i.e., set ^ = id whenever Q 
is inside the scatterer. That is, the relevant phase space is now the unit tangent 
bundle of R d , which will be denoted by T x (M d ). 

The time evolution of an initial particle density / £ L 1 (T 1 (M d )) is 

(3) ft = Lpf 
where is the Liouville operator defined by 

(4) [^/](Q.V):=/(*-*(Q,V)). 

Following Boltzmann's arguments, Lorentz concluded in his 1905 paper that the 
macroscopic time evolution of a particle cloud should, in the limit p —> 0, be 
governed by the linear Boltzmann equation (today also referred to as kinetic Lorentz 
equation), 

° "~^MQ,V)=[ [f t (Q,V Q ) - f t (Q,V)]a(Vo,V)dV , 

Jsl- 1 



(•5) 



8f 



where the collision kernel cr(Vo, V) is the differential cross section of the individual 
scatterer. In the case of elastic scattering at hard spheres we have a(V ,V) = 
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j||Vo — V|| 3_<i . The linear Boltzmann equation describes a random flight pro- 
cess, where a particle moves freely with constant velocity Vq for time t > £ with 
probability e _I/d ^, where Vd is the volume of the (d— l)-dimensional unit ball (the 
total scattering cross section), and is then scattered to velocity V with probability 
<j(V"o, V), again flies with constant velocity V for time t > £ with probability e~ Vd ^ , 
and so on. The crucial observation is that each scattering event is independent of 
the previous one. Thus the process that generates the paths of our random flight 
is Markovian. As we will see, this is different in the case of a periodic scatterer 
configuration. 

The validity of the linear Boltzmann equation was first established rigorously 
for a random, Poisson distributed scatterer configuration by Gallavotti [23]. His 
results were generalized by Spohn |47j to more general random scatterer configu- 
rations and scattering potentials. In 1983, Boldrighini, Bunimovich and Sinai [6 
proved convergence for almost every scatterer configuration drawn from a Poisson 
distribution. As Spohn's work shows, the details of the randomness of the scatterer 
positions is not so essential, and as long as there are no strong correlations, all of 
the above results should remain valid. 

The linear Boltzmann equation has numerous important applications, e.g., to 
neutron transport and radiative transfer, and it is thus natural to ask under which 
circumstances it may fail to provide an accurate description. 



2. The periodic Lorentz gas 

Given a euclidean lattice C C K d of covolume one (e.g., C = Z d ) the periodic 
Lorentz gas is defined as the dynamics of a cloud of non-interacting point particles 
in an array of identical scatterers that are placed at the vertices of the lattice 
C (Fig. [2]). The periodic Lorentz gas has served as a fundamental model in the 
understanding of chaotic diffusion, which emerges in the long-time limit with fixed 
scatterer size. In their pioneering work, Bunimovich and Sinai [8] proved a central 
limit theorem for the long-time dynamics of a particle cloud in two dimensions. 
That is, the long-time evolution converges, in the appropriate scaling limit, to a 
solution of the heat equation. More refined statistical properties that show that 
the dynamics in fact converges to Brownian motion, have recently been established 
in the work of Melbourne and Nichol [41], [42], and the recent paper by Dolgopyat, 
Szasz and Varju [18] . All of the above results assume that the Lorentz gas has finite 
horizon, i.e., there are no infinitely long free flight paths. Without this assumption, 
a central limit theorem still holds, but the diffusion constant is no longer linear 
in time t but diverges as tlogt. This was observed by Bleher [5] and recently 
established rigorously by Szasz and Varju [48]. It is interesting that none of these 
results have so far been extended to dimension d > 3. The arguments given by 
Chernov [TJ] and by Balint and Toth [TJ require non-trivial hypotheses that are 
difficult to establish. 

We now return to the question of the existence of kinetic transport equations 
in the Boltzmann-Grad limit discussed in the previous section. Recall that we are 
interested in the dynamics of a particle cloud 



(6) 



f t (Q, V) = V) = f(-P; l (Q, V)) 
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Figure 2. The Lorentz gas for a periodic scatterer configuration, 
with three distinct particle trajectories. 

in the macroscopic coordinates 0. The estimates by Bourgain, Golse and Wennberg 
0, El] on the distribution of free path lengths already imply that the linear Boltz- 
mann equation does not hold in the periodic set-up; this was pointed out recently 
by Golse [25]. 

The first key result of our joint work with Strombergsson is the proof of the 
existence of a limiting random process for the periodic Lorentz gas [381 Section 1]. 

Theorem 1. Fix a euclidean lattice C. For every t > there exists a linear operator 
L* : L^T 1 ^)) -► L 1 (T 1 (R d )), such that for every f <E L^T 1 ^)) and any set 
A C T 1 (M rf ) with boundary of Lebesgue measure zero, 

(7) lira f [L t p f]{Q,V)dQdV= [ [L* f](Q ,V) dQ dV . 

The operator L* thus describes the macroscopic diffusion of the Lorentz gas in 
the Boltzmann-Grad limit p — > 0. As we shall see however, the family {-L'}t>o does 
not form a semigroup, i.e., 

(8) L s L l f L s+t . 

This is perhaps surprising since {Lp} t >o is indeed a semigroup for every fixed p > 0. 
What is more, in the case of the random scatterer configuration the corresponding 
limiting operators L* also form a semigroup — after all, f t := L* f is a solution of the 
linear Boltzmann equation. The reason for the failure of the semigroup property in 
the periodic setting stems from additional correlations in the lattice, which are lost 
in the macroscopic scaling limit. To keep track of this data, we consider extended 
phase space coordinates (Q, V, £, V+) where (Q, V) £ T 1 (K d ) is the usual position 
and momentum, £ S M+ the flight time until the next collision, and V+ € Sf _1 
the velocity after the next collision. On the microscopic level, the system is now 
over-determined (£ and V + are functions of Q and V), but on the macroscopic 
scale the extra variables are needed. We prove in [351 Section 6] that the particle 
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Figure 3. Two consecutive collisions in the Lorentz gas. 



density ft(Q, V,£, V+) indeed satisfies a generalized linear Boltzmann equation 
(9) 

r\ r\ "I n 

' MQ,V,t,V+)= f t (Q,V ,0,V)po(V 0l V,t,V + )dV . 



dt 



V V 



Q 



The left hand side again corresponds to free transport (note that £ is decreasing lin- 
early with t). The right hand side involves a new collision kernel pq(Vq, V , £, V+), 
given by 



(10) 



Po(V , V, £, V+) = a(V, V+) $ o(£, 6(V, -a(V, V )) 



where cr(y,y + ) is the differential cross section and ^o(£, b(V, V+), — s(V, V"o)j 
the transition probability density to exit with parameter s(V, Vo) and hit the next 
scatterer at time £ with impact parameter b(V,V+); cf. Fig. [3] 

In dimension d = 2 we have the following explicit formula for the transition 
probability [3T)] : 



(11) 
with 

(12) 



6 



T 1 



r 1 



max( to , |z| 



T(x) 



\w + z\ 

' if x < 
a; if < x < 1 
1 if 1 < x, 



The same formula has recently been found independently by Caglioti and Golse 
pTj and by Bykovskii and Ustinov [9], using different methods based on continued 
fractions. 

Our formulas for dimension d > 3 are not as explicit and substantially more in- 
volved, see [37], [3H] for details, and [3D] for asymptotic tail estimates. The formulas 
imply in particular that the collision kernel po (and thus the limiting process) is 
independent of the lattice C on which the scatterers are positioned. Hence any mi- 
croscopic preference for certain directions completely disappears in the Boltzmann- 
Grad limit. 
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The operators 1} in Theorem [T] we were originally interested in can be recovered 
by integrating over the auxiliary variables £ and V+, 



(13) 



where /t(Q, V, £, V_|_) is a solution of the generalized linear Boltzmann equation 
subject to the initial condition 

(14) Jim f t (Q, V, £ V+) - ,g(Q, V)p(V, £, V+) 



wi 



th 



(15) p(V,Z,V+):= [ f a(V ,V)p o (V ,V,e,V + )dV dt;'; 

the latter is a stationary solution of the generalized linear Boltzmann equation. 



3. Why "a generalization" of the linear Boltzmann equation? 

The reason why Q is indeed a generalization of the linear Boltzmann equation 
is the following. As mentioned in the introduction, the linear Boltzmann equation 
corresponds to a random flight process where the time £ until the next collision 
has probability density z^e -1- ' 4 *, where Vd is the volume of the (d — l)-dimcnsional 
unit ball. Furthermore, the probability to exit with parameter s( V, Vo) and hit 
the next scatterer with impact parameter b(V, V+) should be uncorrelated and 
independent of £. We have thus 

(16) sP (£, 6(V, V+), -«(V, V )) = e-""*. 

Substituting in the above the transition density for the random (rather than peri- 
odic) scatterer configuration, we obtain 

(17) Po (Vo, V, £, V+) = a(V, V+) e"^« = p(V, £, V+) 
and 

(18) f t (Q, V, £, V+) = g t (Q, V) a(V, V+) e"^«. 

It is now straightforward to see that ^ yields the classical linear Boltzmann equa- 
tion for gt(Q, V). 



4. Joint distribution of path segments 

The following theorem is the central result in our investigation |38j . It shows that 
the limiting random flight process exists, and that the consecutive path segments 
S\, S%, . . . are generated by a Markov process with memory two. That is, the 
probability for the nth. leg of our path only depends on the previous two path 
segments. 

We will use the notation S := ||S|| _1 S. We set 



(19) B n :={(S u ...,S n )e(R d \{0}) n : S j+1 ? Sj (j = 1, . . . , n - 1)}. 
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Theorem 2. Fix a lattice C and let A be a Borel probability measure on 

which is absolutely continuous with respect to Lebesgue measure. Then, for each 

n G Z>o, and for any set A C M. d x M. nd with boundary of Lebesgue measure zero, 

(20) 

limA({(Q ,V ) eTV~%) : (Q ,Si(Q ,V ; p), ■ • . , S n (Q , V ; p)) G A}) 

= f P^(S 1 ,...,S n )A'(Q Q ,S 1 )dQ dS 1 ---dS n , 

J A 

and where A' is the Radon- Nikodym derivative of A with respect to Lebesgue mea- 
sure. Furthermore, there is a function ^ : i3 3 — > IR>o such that 

n 

(21) PW(S!, ...,S n ) = P< a >(Si, S 2 ) J] *(^_ 2 , Sj-^Sj) 

3=3 

for all n> 3 and all (Si, ... , S„) G 

We in fact also prove a refined version of this theorem, where the initial position 
Q is not random but fixed (on the microscopic scale); see [351 Theorem 1.1]. 
Furthermore, the limiting distributions Pw and W are independent of £ and A. 

We define the probability measure corresponding to (|21[) by 



(22) a4" } (-4):= f P {n) (Si,...,S n )A , (Q ,S 1 )dQ dS 1 ---dS n . 

J A 

Note in particular that (A x K d ) = p^(A). 

5. A LIMITING RANDOM FLIGHT PROCESS 

In Theorem [2] we have identified a Markov process with memory two that de- 
scribes the limiting distribution of billiard paths with random initial data (Q , Vo). 
Let us denote by 

(23) {E(t) : t G K >0 }, 

the continuous-time stochastic process that is obtained by moving with unit speed 
along the random paths Si, S 2 , . . . of the above Markov process with memory two. 
This random flight process is fully specified by the probability 

(24) Pa(S(*i) G Vi, .. .,E(t M ) G T>m) 

that E(t) visits the sets V x , ...,T> M C T 1 ^) at times t = t x ,.. .,t M , with M 
arbitrarily large. To give a precise definition of (23 1 set To := 0, T n := Y^j=i ll^ill; 
and define the probability that S(£) is in the set T>\ at time t\ after exactly n\ hits, 
in the set T> 2 at time t 2 after exactly n 2 hits, etc., by 



(25) 



G 2?i, . . . , S(£m) G 2? Af, T ni < ti < T ni+ i, . . . , T„ M < tu < T nM+ i) 
:= Ma" +1) • • • , S„+i) : S nj . fo) G 23 i5 T nj < t s < T n . +1 (j = 1, . . . , M)}) 

with n := (raj., . . . , tt,m), n := max(ni, . . . , um), and 

(26) E n (t) := S rf + (t - T n )S n+1 ,S n+1 
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The formal definition of (|23l) is thus 



(27) V A (E(t 1 )eV 1 ,...,E(t M )€V 



M 



]T Pi n) (s(ti)GPi,...,s(( M )ei) j; 



and T ni <ti< T 7ll+ i, . . . ,T nu < tu < T nM +i). 
The following theorem shows that the Lorentz process (l2| converges to the sto- 



chastic process ( 23 1 as p — > 



Theorem 3. Fix a lattice C and let A be a Borel probability measure on T 1 (R d ) 
which is absolutely continuous with respect to Lebesgue measure. Then, for any 
t\, . . . , tM S K>0; and any subsets T>\, . . . , T>m C T 1 (M d ) with boundary of Lebesgue 



measure zero. 



(28) limA({(Q ,y )eT 1 (p d - 1 /C p ): 

{Q{tx),V{tx)) eV 1 ,...,{Q{t M ),V(t M )) eV M }) 

= P A (H(tO e Vu..., E{t M ) GPjf). 
T/ie convergence is uniform for t±, . . . , 4m * n compact subsets ofR> . 

Theorem [3] follows from Theorem [2] the main ingredient in the proof is an esti- 
mate that shows that it is unlikely to have many collisions in any fixed time interval 
[38j Section 5]. 

The generalized linear Boltzmann equation Q can now be interpreted as the 



Fokker-Planck-Kolmogorov equation of the stochastic process ( 23 1 , and its valid- 
ity follows from Theorem [3] by standard arguments from the theory of stochastic 
processes |38l Section 6]. 



6. The distribution of free path lengths 

To give a detailed account of the proof of Theorem [2] would go beyond the scope 
of this lecture. I however hope to be able to explain the key idea, which is already 
clearly visible in the case n = 1, i.e., the distribution of free path lengths. 

In order to explain the proof, it will be more convenient to return to the original 
microscopic coordinates (q,v), where the free path length diverges at the rate 
p-Cd-i). 

Let us denote by JC p C K d the complement of the union of all scatterers. The 
free path length for the initial condition (q,v) € T 1 (/C p ) is defined as 

(29) n(q,v;p)=mf{t>0:q + tv<£lC p }. 

That is, Ti(q, v; p) is the first time at which a particle with initial data (q, v) hits 
a scatterer. We also include the (somewhat artificial) case when q £ £, i.e., when 
the particle starts at the center of a scatterer; in this case we think of K. p C K d as 
the domain obtained by removing all scatterers except the one centered at q. The 
following is the main result of [37j . 

Theorem 4. Fix a lattice C of covolume one, let q E R d , and let X be a Borel 
probability measure on Sf^ 1 absolutely continuous with respect to Lebesgue measure. 
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Then there exists a C 1 function Fc, q on R>o such that, for every £ > 0, 

(30) lim A({« S S?- 1 : ^nfe v; p) < £}) = F c<q (0- 

The distribution of the free path lengths in the Lorentz gas was already investi- 
gated by Polya, who rephrased the problem in terms of the visibility in a (random 
and periodic) forest (43]. The problem of the limiting distribution in dimension 
d = 2 was recently solved by Boca and Zaharescu [5J in dimension d = 2 in the 
case when q is either random or located at a lattice point; see also their earlier 
work with Gologan [3j , and the paper by Calglioti and Golse [TO] . Previous work in 
higher dimension d > 3 includes the papers by Bourgain, Golse and Wennberg [7], 
|24j who provide upper and lower bounds on the tail of the distribution of free path 
lengths. More details on the existing literature can be found in the survey |25j. 

7. The space of lattices 

The new idea in the joint work with Strombergsson [37] is to translate the prob- 
lem of the free path length into a question about the dynamics on the space of 
lattices. The advantage of this approach over previous attempts is that the tech- 
nique extends naturally to arbitrary dimension, and that the limiting distributions 
have a canonical interpretation as the distribution function of random variables on 
a beautiful geometric object. 

A euclidean lattice £ C ff 1 of covolume one can be written as C = Z d M for 
some M € SL(rf, R), where Z d is the standard cubic lattice. Since SL(d, Z) leaves 
Z d invariant, the homogeneous space X\ — SL(d, Z)\ SL(d, R) parametrizes the 
space of lattices of covolume one. Similarly, let ASL(rf, K) = SL(d,K) k R d be the 
scmidirect product group with multiplication law 

(31) (M, £)(M', O = (MM 1 , £M' + £')• 
An action of ASL(d, R) on R d can be defined as 

(32) y^y(M,S):=yM + S. 

Each affine lattice (i.e. translate of a lattice) of covolume one in R d can then be 
expressed as 1 d g for some g € ASL(d, R), and the space of affine lattices is then 
represented by X = ASL(d,Z)\ ASL(d,R) where ASL(d,Z) = SL(d,Z) x 1 d . We 
denote by \x\ and \i the Haar measure on SL(d, R) and ASL(d, R), respectively, 
normalized in such a way that they represent probability measures on X\ and X. 
We are interested in the A-measure of velocities with free path lengths at most 

(33) \({v eSt 1 : p^rxiq^p) <£}). 

This is approximately the same as the A-measure of directions such that a cylinder 
q+Z(v : p - ^ -1 ^! p) of length p _ ( d_1 )£ and radius p, pointing in direction v contains 
at least one lattice point (cf . Fig. [1]) : 

(34) « \{{v € St 1 : Z d M Dq + Z(v, p'^k, p) + «})■ 

The approximation comes from the fact that our cylinder should have spherical 
caps of radius p on each end; it is easy to show however that the A-measure of 
v that have a lattice point in these caps is vanishingly small, as p — > (see [37] 



10 



JENS MARKLOF 



OOOOOOOOOOO 

ooooooooooo 
000,00000000 
ooooooooooo 
ooooVooooooo 
\ cM^ 0000 
000 ^oj'o OOOO 
0000 o\c5^o 0000 

00000 \pl 00000 
A^ 

ooooodooooo 
000000 \ o OOOO 

OOOOOOOOOOO 
OOOOOO o\ OOOO 
OOOOOOOOOOO 

00000000.000 
OOOOOOOO \o 00 
ooooooooooo 
ooooooooooo 
ooooooooooo 




'if 






d- 










: * .* * •'.* • 'frtf-vi: -V • 


' ' v.v. 



Figure 4. Renormalization of the periodic Lorentz gas. Step 1: 
Replace the question on the number of intersections of a line with 
scatterers by a lattice point counting problem in a long stretched 
cylinder. Step 2: Rotate cylinder and lattice so the cylinder lies 
horizontally. Step 3: Apply a diagonal linear transformation that 
maps the cylinder to the p-independent reference cylinder 2(£). 

Section 4.1] for details). The next step is to shift lattice and cylinder by — q, and 
then rotate by K (v) G SO(d) such that K(v)v — ei, where e% = (1,0, . . . , 0): 

(35) = \({v G St 1 : {1 d M - q)K{v) n 2(e 1 ,p^ 1 ^, p) + 0}). 

The cylinder now lies in the ei-direction. We apply the linear transformation 
D(p) = diag(p d_1 , . . . ,p _1 ) which transforms the long and thin cylinder into 
better proportions: 

(36) = \({v G St 1 ■ {% d M - q)K{v)D{p) n z(0 + 0}), 
where 

(37) Z{i) := Z(ei,f, 1) = • ■ ■ , x d ) G R d : < x x < £, \\(x 2 , . . . , x d ) \\ < l}. 

Although the above linear transformations seem trivial, we have achieved a different 
perspective on the problem: Rather than counting lattice points in a long thin 
cylinder (which looks hard) we now count in a well proportioned object. The 
prize we have paid is that our original lattice C = Z d M has changed to (Z d M — 
q)K(v)D(p). Hence we are moving through the space of (affine) lattices, as p — > 0, 
and may now employ ergodic theoretic methods to understand the averages over v 
with respect to A. 
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8. EQUIDISTRIBUTION IN THE SPACE OF LATTICES 

We begin with the simplest case, q £ C, i.e., without loss of generality q = 0. 
The right translation 

(38) X X ^X X , SL(d,Z)M^ SL(d,Z)M$* 



by the element 

(39) $* 



e -(d-i)t o 
*0 e 4 !^! 



defines a flow on the homogeneous space X\ = SL(d, Z)\ SL(d, K). This flow has 
many good chaotic features: it is ergodic, mixing and partially hyperbolic. 
The horospherical subgroups generated by 

(40) n + (x)= J), «_(*)= J 

respectively, generate the stable and unstable horospherical subgroups of the flow. 
Using the mixing property, a standard argument (see e.g. |2T], [33]) shows that aver- 
ages over the unstable horosphere become asymptotically (t — > oo) equidistributed 
in Xi with respect to 

Theorem 5. Let X be a Borel probability measure on absolutely continuous 

with respect to Lebesgue measure. Let f : X x — ► K be bounded continuous, M £ 
SL(d,R). Then 

(41) lim / /(Mn_(x)$*)dA(») = / f(M')dm{M'). 



t- 



Alternative proofs of this theorem can be obtained by using harmonic analysis 
(which is particularly feasible when M = 1, i.e., the horosphere is closed) or Ratner's 
theory, which we will revisit below. 

In view of ( 36 1 we are interested in the distribution of the orbit 

(42) SL(d 1 Z)\{SL(d,Z)MK(v)^ t : v £ S^ 1 } 

in X\ , as t — » oo (set t = log 1/ p). By using the fact that this orbit is close to an 
unstable horosphere, we can exploit Theorem [5] to deduce equidistribution also in 
this case; we refer the reader to [37] Section 5] for details. 

Theorem 6. Let X be a Borel probability measure on S^ _1 absolutely continuous 
with respect to Lebesgue measure. Let f : X\ — > K be bounded continuous, M £ 
SL(d,M). Then 

(43) lim / f(MK{v)&)d\{y) = f f(M')dfj,{M'). 

t^ocJ s d-l J Xi 

This theorem thus states, that in the limit t — > oo we can replace the A-average 
over v by an average of the entire space of lattices. This yields (modulo some 
technicalities) the proof of the limit law for the free path length, Theorem |1J in the 
case q £ C, plus a formula for the limit distribution: 

(44) F Ct0 (0 = fii({M £ X 1 : Z d M n ^ 0}). 

That is, the limit distribution of the free path length for a particle emerging from the 
center of a scatterer is given by the probability that a random lattice intersects the 
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cylinder Z '(£) m a£ /easf one point. Note that the limit distribution -Fb(£) := i*£,o(£) 
is independent of the lattice C and of A. 

Instead of particles emerging from a lattice point we can also consider initial 
conditions on the boundary of a scatterer; this leads to different limit distributions 
and is one of the crucial steps in the proof of the existence of the limiting random 
flight process described in Theorems [2] and [3] 

Let us now turn to the case q £ QC. The right translation 

(45) X^X, ASL(d,Z)g^ ASL(d,Z) 5 $ t 

by the element 

,t /Ye-**- 1 )* 



(46) *=IA * e«l<_J' 

now defines a flow on the homogeneous space X = ASL(d, Z)\ ASL(d, R). In anal- 
ogy with the above, we set 

(47) n-(x)=(£ d - 1 Y)>° 

The crucial difference is now that n_(x) no longer generates the full unstable horo- 
sphere for the flow <!>', and hence the mixing argument is no longer sufficient. 
Instead, we need to employ Ratner's classification of measures that are invariant 
under unipotent actions, which in the present case is given by the right action of 
n-(x). We can in particular exploit a very useful theorem of Shah (46] to show 
the following. (See [37J Section 5] for details, and [3B] for a general introduction to 
applications of Ratner's theory to problems of this kind.) 

Theorem 7. Let A be a Borel probability measure on absolutely continuous 

with respect to Lebesgue measure. Let f : X —> R be bounded continuous, a G 
R d \Q d , Me SL(d,R). Then 

(48) lim / f((l,a)(M,0)n-(x)&)d\(x)= [ f(g)dn(g). 



t — >oo 



X 



Using the same approximation argument that led to Theorem [6] we deduce: 

Theorem 8. Let X be a Borel probability measure on S^ _1 absolutely continuous 
with respect to Lebesgue measure. Let f : X — > K be bounded continuous, a G 
R d \ Q d , M € SL(d, R) . Then 

(49) lim / f((l, a )(M,0)(K(v),0)¥)d\(v)= [ f(g)d^g). 

This proves Theorem [4] in the case q £ QC, and yields the formula 

(50) F c , q (0 = x) € X : (Z d M + x) n Z(0 ? 0}). 

Hence, the limit distribution of the free path length for a particle starting at a generic 
position is given by the probability that a random affine lattice intersects the cylinder 
Z{£) in at least one point. Here, F(£) := Fc, q (^) is evidently independent of C, q 
and A. 

I omit the discussion of the remaining q £ QC — the arguments are analogous to 
the above [571. 



KINETIC TRANSPORT IN CRYSTALS 



13 



Only in dimension d = 2 we have been able to turn the formulas ( 44 1 and ( 50 ) 



for -Fb(£) and F(£) into explicit functions of £. The easiest approach is to integrate 
our formula for the transition probability (11) (see [39] for details), 



(51) 



(52) F(£) = / / / / M£,",w,z)dwdzd£,"d£,'. 

JO J£' J-l J— I 

The resulting explicit expressions coincide with those obtained, using different 
methods, by Boca, Gologan and Zaharescu [4], and by Boca and Zaharescu [5], 
respectively. 



9. ASYMPTOTICS 

The geometry of the space of lattices is significantly more complicated in di- 
mension d > 3, and it seems extremely hard to obtain any explicit formulas. It 
is possible, however, to describe the asymptotic tails of our limit distributions, by 
observing that when £ > is very large (or small), then the lattices that contribute 
to -Fb(£) and F(£) must have at least one very short basis vector. That is, all 
integration is restricted to the cusps of the spaces X\ and X, respectively, whose 
geometry is simpler. Using tools from the geometry of numbers we can show that 

m 

(53) Fo(Q — 1 f° r £ sufficiently large, 
and 

(54) ^o(£) = ^y£ + 0(£ 2 ), £-0, 

where = Y((d+i)/2) ^ s ^ ne v °l ume of the (d— l)-dimensional unit ball. Similarly, 



d-l 

7T 2 



(55 » m = 1 -¥mmw> rl+oir '~ ,) ' ( ~°°- 

and 

(56) F(£) = ^£ + 0(£ 2 ), £^0. 

We also obtain asymptotic formulas for the collision kernel <?o(£, w, z) in the limits 
£ — > and £ — > oo, see [40 for details. 



10. Outlook 

The techniques outlined above do not necessarily require that the scatterers are 
rigid spheres. It is sufficient to assume that the scattering map is dispersive; a 
Muffin-Tin Coulomb potential would be a good example [3H] . A key hypothesis of 
our approach is however is that the interaction region of the scattering process is 
finite, so that the test particle moves along straight lines for most of the time. This 
assumption is no longer valid in the case of long-range potentials. Provided the 
renormalization approach can be modified accordingly, it seems feasible to gener- 
alize our studies to crystals with long-range potential, at least for sufficiently fast 
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power-like decay. A result by Desvillettes and Pulvirenti [15] achieves this objec- 
tive (with some additional technical assumptions) in the case of random scatterer 
configurations; cf. also the work by Poupaud and Vasseur [33]. 

A related problem is to consider different scaling limits for compact potentials, 
where the strength of the potential is reduced, and at the same time the scatterer 
density rescaled suitably to achieve a non-trivial limit. In this case grazing collisions 
become important, and one expects a different kinetic equation for the macroscopic 
dynamics; cf. Desvillettes and Ricci 1 GJ for the corresponding result for a random 
scatterer configuration — here the limiting kinetic equation is the classical Fokker- 
Planck equation. 

The renormalization approach we have developed for the Lorentz gas assumes 
that the scattering process is the same for each scatterer. It is possible to remove 
this assumption, as long as the resulting (global) potential of the crystal is still 
periodic, or in the case of quasicrystals. A natural question leading on from this 
is whether one can describe the kinetic equations in the case of a (quasi-)crystal 
with random defects, i.e., a scatterer is removed from a lattice site with probability 
p, with < p < 1. Such a set-up will lead to an interesting variation of the 
renormalization method, since the modular invariance of the crystal will be replaced 
by a modular invariance in distribution (we assume p is the same for each scatterer) . 
It is no surprise that the limiting case p — > 1 leads back to the linear Boltzmann 
equation |45j . 

A further important case that has been extensively studied for stochastic Lorentz 
gases is the dynamics in the presence of electro-magnetic fields [17] . In the case of 
constant magnetic or electric fields, the above problem has some connection with 
beautiful, basic number-theoretic questions on the distribution of lattice points near 
circles or parabolas, respectively. 

Finally, an obvious challenge is to adapt our renormalization approach to the 
quantum mechanical problem, and derive a corresponding quantum kinetic equa- 
tion, in suitable weak-coupling or low-density limits discussed for random scatterer 
configurations by Erdos and Yau [2DJ , and Eng and Erdos [ID] ■ It should be stressed 
that of course the quantum theory of electrons in fixed periodic potentials is well 
understood, and that the lattice symmetry allows for a wealth of techniques in the 
spectral analysis (Floquet-Bloch decomposition). In particular the KKR method 
is an extremely useful tool, and has been applied to the semiclassical analysis of 
the symmetry-reduced periodic Lorentz gas (the Sinai billiard) in connection with 
quantum chaos, see Berry [2], and (for the case of small scatterers) Dahlqvist and 
Vattay [H] . 

The central idea of the work presented here is to exploit the dynamics of flows 
on the space of lattices. This technique has proved to be extremely powerful also 
in other applications in mathematical physics, including KAM theory [27], [25], 
arithmetic quantum unique ergodicity [3DJ and the Berry- Tabor conjecture [22], 
[32], [33], [35]. Flows on the space of lattices (and more general homogeneous 
spaces) should be viewed as the higher-dimensional generalization of the classical 
continued fraction algorithms — I expect more striking applications in the future. 
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